Integrable models for asymmetric Fermi superfluids: Emergence of a new exotic pairing phase 
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We introduce an exactly-solvable model to study the competition between the Larkin-Ovchinnikov-Fulde- 
Ferrell (LOFF) and breached-pair superfluid in strongly interacting ultracold asymmetric Fermi gases. One can 
thus investigate homogeneous and inhomogeneous states on an equal footing and establish the quantum phase 
diagram. For certain values of the filling and the interaction strength, the model exhibits a new stable exotic 
pairing phase which combines an inhomogeneous state with an interior gap to pair-excitations. It is proven that 
this phase is the exact ground state in the strong coupling limit, while numerical examples in finite lattices show 
that also at finite interaction strength it can have lower energy than the breached-pair or LOFF states. 
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Superconducting and Fermi superfluid phenomena have 
been a subject of fascination since their discovery. Both phe- 
nomena are direct macroscopic-scale manifestations of quan- 
tum physics, with the electric charge of their relevant mi- 
croscopic constituents being the crucial factor differentiating 
them. Interest in their various fundamental aspects has in- 
creased recently because advances in the field of ultracold 
atomic Fermi gases QJ 0] are leading to new experimental 
probes to investigate unexplored territory, with consequences 
in condensed matter as well as high energy physics (e.g., the 
physics in the core of neutron stars). 

Of particular importance is the nature of their ground states 
(GSs) under various external conditions since novel thermo- 
dynamic phases might show up. The present manuscript stud- 
ies the relative vacuum stability of a two-species fermion gas 
as a function of the pairing interaction strength g and differ- 
ent species population. Without loss of generality we denote 
the two species as a and b with densities p a ^)- Differences 
among the two species could be related to their masses m Q (6), 
spin or hyperfine states. Asymmetry, in general, makes pair- 
ing less favorable and questions about the nature of the result- 
ing competing phases might arise. To this end, we introduce a 
model system that will prove to be exactly solvable and which 
displays the competition between the Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) @|, breached-pair (or Sarma) 001- de- 
formed Fermi-surface superfluid @|, and segregated phases 
0]. Its quantum phase diagram as a function of the asymme- 
try density Sp = (p a —pb)/ (Pa+Pb) and the coupling strength 
g has been recently studied, at the mean-field level, in the one- 
channel 1 8] and two-channel models |;9L In our exact solution, 
albeit in a finite-size lattice, we find different regimes of sta- 
bility for these various phases. A key result is the prediction of 
a new exotic inhomogeneous phase characterized by a partic- 
ular center-of-mass momentum of the condensed pairs. This 
phase is the exact GS in the large-g limit. <=i V 

Consider N a and Nf, fermionic atoms confined to a D- 
dimensional box of volume V, i.e. p a ^ = N a n,\/V, with where \v T ) = Y[ 



ents to study the competition between the various phases 
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where a k (& k ) creates a particle of type a(b) with momentum k 
and n£ = a k a k , n k = b\p k . The pairing interaction scatters 
pairs with center-of-mass momentum Q and "band" energies 
eg = ek/2m Q (a = a,b), with ek representing an arbitrary 
dispersion (including a non-rotational-invariant one 01). 

The quantum integrability and exact solvability of the 
Hamiltonian ([0 can be derived using an su(2) algebra 
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and a second, independent, realization of su(2) 
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These two mutually commuting algebras are often referred to 
as charge and spin su(2) realizations, respectively. Usin g th e 
algebraic techniques of the Richardson-Gaudin model jlQIl . 
one can write down a complete set of integrals of motion 
i?T Q , W ith [Rl Q ,RZ tQ ] =0 (for T,T' = r,S): i£ >Q = 
T k,q + 2 5TEk'(^k) X kk' T k ,q • f k ',Q, where X£ k , = 
V(Wk Q) ~~ ^(k' q))' W ^ m arbitrary functions rj T depending 
upon k and Q, and gj are the coupling constants. Their com- 
plete set of eigenvectors are of the form 
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periodic boundary conditions and g < 0. (The exact solvabil- 
ity of the problem is not restricted to these latter conditions 
but for notational convenience we will specialize to this case). 
The following model Hamiltonian contains the right ingredi- 
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is a quasispin vacuum state de- 



fined by T^q \v J ) = 0, and T£ q \v 1 ) = d k q |^ t ), with 
dJ.Q = (2^k q — O k ,Q)/4, rik,Q = 2, and the seniority quan- 
tum number v^q = 1,0, which for the su(2) pair algebra Q 



counts the number of unpaired fermions. The complex spec- 
tral parameters Ej satisfy the set of non-linear equations 
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To simplify matters, and because our goal is to show ex- 
act solvability of Q, we will only consider dynamics in 
the charge space (i.e., 175 = 0, g T = g, dropping the la- 
bel T). The total number of atoms is N = N a + Nb = 
2M + v, where M is the number of atom pairs and v the 
number of unpaired ones. Consider now the linear com- 
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C, where C 
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g (N — L) /2 — g (N — L). Comparing H T with Q we im- 
mediately see that they differ in the kinetic term. Making 
use of the spin su(2) algebra by adding a term of the form 
2 Ek £(k,Q)>Sk Q' it leads to (up to an irrelevant constant) 

H = yi(V(k,Q) + £(k,Q))K+Q + (V(k,Q) -£(k,Q))™_ k ] 
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Identifying r ?(k , Q) = §[£k +Q - 

e^ k ], we get the Hamiltonian of Eq. Q, after constraining the 
vectors k + Q and k to be in the same set. The eigenvalue E 
corresponding to the solutions of Eqs. (jSJl are given by 
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where denotes the number of unpaired a particles in the 
state with momentum k. The space dimensionality of the 
problem enters through the band dispersion eg, and the ef- 
fective degeneracies d k ,Q in the exact solution The latter 
are in turn defined by ?7(k,QV Assuming space-inversion sym- 
metry (eg = e™ k ), the degeneracies Slk.Q count the number 
of states [k, Q] with the same value of rjtk q). 

Our exactly-solvable model is valid for arbitrary Q values. 
The Q = limit restores the homogeneous BCS phase giving 
rise to a breached-pair phase in terms of 0-momentum pairs, 
while a finite value of Q gives rise to the LOFF phase with 
Q-momentum pairs. For an asymmetric system with an ex- 
cess of the a species (N a > Nb), the atoms fill the lowest 
states up to k^ with |k^| > |k^| at weak coupling. When 
the interaction is switched on, several possible states com- 
pete to determine the absolute GS. The position of the un- 
paired atoms, defining the seniority quantum numbers z^k,Q, 
block the available states from scattering pairs of atoms effec- 
tively reducing the degeneracies to d k ,Q- When Q = the 
equations reduce to the well-known Richardson model with 
blocked states fiol [Till . In general, configurations are iden- 
tified by their g — > limit, with specific pair and seniority 
occupations. They can be categorized as follows: 

• asymmetric BCS (aBCS): Q = 0, a and b particles fill 
their lowest orbitals up to their corresponding Fermi levels. 
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FIG. 1: Three Q = configurations, asymmetric BCS (aBCS), 
breached A, and breached B, as described in the text for a D — 1 
lattice with 12 sites, N a = 9 and Nb = 5 atoms. The left vertical 
axis displays the single-particle energies e£ — — 2cosk (a = a, b) 
while the right one shows the corresponding momenta k. An allowed 
pair-scattering process is indicated in each case with arrows. 
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FIG. 2: Two possible Q = ir/3 configurations. The first one (LOFF) 
corresponds to a Fermi sea. A possible pair-scattering process is 
indicated with arrows. We show a breached LOFF configuration as a 
second example. The k momentum values are the same as in LOFF. 



• breached A: same as aBCS, but the unpaired a particles 
move up in energy such that pairing correlations can develop 
around k£, . 

• breached B: same as aBCS, but the unpaired a particles 
move down in energy such that pairing correlations can de- 
velop around k p . 

• LOFF: finite Q, a and b particles fill their lowest orbitals up 
to their corresponding Fermi levels. 

• breached LOFF: finite Q, but now some of the unpaired a 
particles move to allow more pairing correlations. 

We illustrate some of these states by using a D = 1 lattice 
with L = 12 modes as an example. In Fig. ^ we show the 
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level scheme for a system with N a = 9 and N b = 5. In the 
first column, labeled aBCS, the excess of a atoms occupy the 
states between k^, and k F completely blocking these states. 
The corresponding b states are represented by a dash line. Pair 
scattering can only occur between states below k^ and states 
above k F , as indicated in the figure. It is worth noting that 
the excess of a atoms could be located in any configuration in 
k-space. In the second column we display a breached-pair 
superfluid state ^ 0] (Breached A) where the unpaired 
a atoms are promoted to higher-energy states to leave some 
space around k b F for pairing. A possible pair-scattering pro- 
cess is indicated in the figure. Alternatively, the blocked states 
could be moved down for the pair scattering to take place 
around k. F as shown in the third column of the same figure 
(Breached B). The relative stability of each one of these pos- 
sible states will depend upon the competition between the ki- 
netic energy and the pairing interaction. 

Figure[2]displays two examples of a LOFF state configura- 
tion. Here we assume a momentum Q that exactly matches 
the two Fermi energies (Q = k F — k^, = 7r/3). In the first 
example, corresponding to the first two columns of Fig. |2] the 
atoms occupy the lowest single-particle energies e£ defining 
a configuration which is expected to be the lowest LOFF state 
at weak coupling. The numbers within the circles indicate the 
Q = 7r/3 momentum pairs. The unpaired a atoms, displayed 
with a black circle, block the corresponding states of the b 
atoms. A possible breached LOFF configuration is shown as 
a second example. 

We will now explore the competition between the possi- 
ble phases in a numerical example for D = 2. We assume a 
square lattice with dispersion eg = — 2(cos(k x ) + cos(k y )), 
with units chosen such that k x and k y are multiples of 2-k/L, 
where L is the linear size of the lattice. Being the disper- 
sion equal for both atomic species, we are excluding an asym- 
metry in the masses or a deformed Fermi surface. Prelimi- 
nary results for asymmetric masses do not show qualitative 
differences with the results presented below, except that for 
m a > m b breached A will be more stable than breached B 
because particles of type a will require less kinetic energy to 
shift to higher momenta than particles of type b. 

We discuss first the limiting cases of a very weak or very 
strong pairing interaction. If the interaction strength g is much 
smaller than the level spacing, then the full problem reduces 
to a pairing problem for each level separately, with the cou- 
pling between levels entering at order g 2 . The leading order 
is the single-particle energy. This means that the GS fills the 
lowest single-particle orbitals up to the Fermi levels of species 
a and b, respectively. This leaves no room for the breached- 
pair phase in the GS. aBCS and LOFF are degenerate to lead- 
ing order. This degeneracy is lifted at first or second order in 
g. The pairing interaction favors open shells, hence it prefers 
a non-zero value of Q when the valence shell is completely 
filled, while it might prefer Q = for open valence shells. 
For the 6x6 model studied here, we found that aBCS dom- 
inates the weak limit for N a — N b < 4 at quarter filling and 
for N a - N b < 10 at half filling. 

In the strong-coupling limit one can expand Eq. (JSJl in terms 
of g~ x . In this way, the asymptotic GS of Q can be analyti- 
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FIG. 3: Quantum phase diagram for a 6 x 6 lattice at half filling; the 
inset displays the quarter-filling case (aBCS: yellow, LOFF: green, 
breached B: red, breached LOFF: blue, dark blue for Q = (n, tt)). 
The shaded area indicates the transition from the normal to the su- 
perfluid phase llol . 

cally determined fl2ll . The resulting GS energy is given by 

E = -2gN b (2fl+l-N b )+J2^>i 

k 

+ 2N b ^ Q ^ +0(g-^ (8) 

with n = ^k.Q- Upon inspection one finds that for the 
lattice model considered here, the lowest possible value for Eq 
will occur when Q = (tt, tt). In that case all ?7(k,Q) vanish, 
and the first line of Eq. (jSJi becomes an exact expression with 
the excess a atoms occupying the lowest single-particle states. 
One can describe this regime as an extreme breached LOFF 
state. This result is exact in the limit of |(/| much larger than 
the bandwidth, which is an unphysical assumption. However, 
it indicates that at some finite value of g a transition to an ex- 
otic inhomogeneous phase must occur, combining a breached 
configuration with a non-zero value for Q. We find such con- 
figurations to have a lower energy than the aBCS, breached A 
or B or LOFF configurations at interaction strengths as weak 
as g = —0.1, which might be realizable in a physical setting. 

We have studied model Q numerically on a 6 x 6 lattice, 
with 18 (quarter-filled) or 36 (half-filled) particles, distributed 
over a and b states. Finding the optimal configuration for the 
unpaired particles turned out to be highly non-trivial because 
of the large number of possibilities. We addressed this prob- 
lem using a quantum Monte Carlo technique fl3ll that pro- 
vided a number of candidate GS configurations for various 
values of g and Q. Starting from the g = configurations, the 
solution of Eq. (0 was then obtained by slowly increasing the 
value of q, and by applying the iteration techniques explained 
in Ref. 11411 . In this way we evaluated the exact energies for 
each configuration up to g = —0.5, and we were able to de- 
termine the GS and the exact transition points. One can see 
in Fig. [3]that exotic configurations such as LOFF or breached 
LOFF can have a lower energy than the aBCS state. There is 
a subtle competition between LOFF and the various breached 
BCS states, and both phenomena appear simultaneously in the 
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FIG. 4: Occupation numbers in momentum space for various con- 
figurations with particle numbers N a = 22 and Nt = 14, at 
g = -0.5, for a 6 x 6 lattice (LOFF:Q = (tt/3,0); breached 
LOFF:Q = (tv/3, 7t/3)). Occupation numbers for the unpaired par- 
ticles have been symmetrized over all possible orientations. 



emergent breached-LOFF regime that dominates the phase di- 
agram at larger asymmetries and interaction strengths. The 
recently observed 11511 normal state region at weak coupling 
could be qualitatively determined using the techniques dis- 
cussed in lUal . Our approach also allows computation of the 
occupation numbers in momentum space. They are derived 
from the integrals of motion using the Hellman-Feynman the- 
orem. Figure|4]shows results for a selected number of config- 
urations, corresponding to the lowest-lying states of the half- 
filled model at N a - N b = 8, to illustrate the aBCS, breached 
A and B, LOFF and breached LOFF phases at g = —0.5. 

In summary, we presented a one-channel exactly-solvable 
model that admits several homogeneous and inhomogeneous 
phases depending upon the relative strength between kinetic 
and pairing interactions, and the difference in the number of 
atomic species (given a fixed total number of atoms). The 
inhomogeneous phases (LOFF) show up as soon as the differ- 
ence in Fermi momentum between the two species becomes 
commensurate with the unit lattice momentum. A most sig- 
nificant result is the prediction of a new exotic phase which 
combines pairs with definite momentum and breached super- 
fluidity/superconductivity, that we dubbed breached LOFF. 
We expect this new phase to be the stable ground state at 
large interaction strengths for fixed, asymmetric particle num- 
bers. These phases can be experimentally differentiated in 
time-of-flight measurements of the molecular velocity, after 
sweeping the system through the BCS-to-BEC crossover re- 
gion (ll 01 . The momentum distribution of unpaired fermions 
may distinguish the various exotic phases discussed here. The 
present analysis can also be extended to a two-channel inte- 
grable model with the explicit treatment of the Feshbach res- 
onance 1 9] in a similar way as in Ref. fl7ll . 

We acknowledge discussions with J. Carlson, W. V. Liu, 
C. Lobo, S. Reddy and D. Van Neck. JD, SR and KVH 
acknowledge financial support from the Spanish DGI under 
grant BFM2003-05316-C02-02, and the FWO - Flanders. 



[1] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 92, 
040403 (2004); M. W. Zwierlein et. al, Phys. Rev. Lett. 92, 
120403 (2004). 

[2] G. Ortiz and J. Dukelsky, Phys. Rev. A 72, 043611 (2005). 
[3] A. I. Larkin and Yu. N. Ovchinnikov, Sov. Phys. JETP 20, 762 

(1965); P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964). 
[4] G. Sarma, J. Phys. Chem. Solids 24, 1029 (1963). 
[5] W. V. Liu and F. Wilczek, Phys. Rev. Lett. 90, 047002 (2003). 
[6] H. Muther, and A. Sedrakian, Phys. Rev. Lett. 88, 252503 

(2002). 

[7] P. Bedaque et al, Phys. Rev. Lett. 91, 247002 (2003). 

[8] A. Sedrakian et al, Phys. Rev. A 72 013613(2005); C. -H. Pao 

et al, Phy s. Rev. B 73, 132506 (2006); D. T. Son, and M. A . 

Stephanov, cond-mat/0507586 K. Yang, cond-mat/0508484 
[9] D. E. Sheehy, and L. Radzihovsky, Phys. Rev. Lett. 96, 060401 

(2006). 



[10] J. Dukelsky et al, Rev. Mod. Phys. 76, 643 (2004); G. Ortiz et 

al, Nucl. Phys. B 707, 421 (2005). 
[11] R. W. Richardson, Phys. Lett. 3, 277 (1963); Nucl. Phys. 52, 

221 (1964). 

[12] R. W. Richardson, J. Math. Phys. 18, 1802 (1977); E.A. 
Yuzbashyan, A.A. Baytin, and B.L. Altshuler, Phys. Rev. B 68, 
214509 (2003). 

[13] S. M. A. Rombouts, K. Van Houcke, and L. Pollet, Phys. Rev. 

Lett. 96 180603 (2006). 
[14] S. Rombouts, D. Van Neck, and J. Dukelsky, Phys. Rev. C 69, 

061303 (2004). 

[15] M. W. Zwierlein etal, Science 311, 492 (2006); G. B. Partridge 

etal, ibid. 311, 503 (2006). 
[16] J. Dukelsky et al, Phys. Rev. Lett. 88, 062501 (2002). 
[17] J. Dukelsky etal, Phys. Rev. Lett. 93, 050403 (2004). 



4 



